Molecular Mechanisms Underlying Medium-Chain Free Fatty Acid-Regulated Activity of the Phospholipase PlaF from Pseudomonas aeruginosa.

PlaF is a membrane-bound phospholipase A1 from Pseudomonas aeruginosa that is involved in remodeling membrane glycerophospholipids (GPLs) and modulating virulence-associated signaling and metabolic pathways. Previously, we identified the role of medium-chain free fatty acids (FFAs) in inhibiting PlaF activity and promoting homodimerization, yet the underlying molecular mechanism remained elusive. Here, we used unbiased and biased molecular dynamics simulations and free energy computations to assess how PlaF interacts with FFAs localized in the water milieu surrounding the bilayer or within the bilayer and how these interactions regulate PlaF activity. Medium-chain FFAs localized in the upper bilayer leaflet can stabilize inactive dimeric PlaF, likely through interactions with charged surface residues, as has been experimentally validated. Potential of mean force (PMF) computations indicate that membrane-bound FFAs may facilitate the activation of monomeric PlaF by lowering the activation barrier for changing into a tilted, active configuration. We estimated that the coupled equilibria of PlaF monomerization-dimerization and tilting at the physiological concentration of PlaF lead to the majority of PlaF forming inactive dimers when in a cell membrane loaded with decanoic acid (C10). This is in agreement with a suggested in vivo product feedback loop and gas chromatography-mass spectrometry profiling results, indicating that PlaF catalyzes the release of C10 from P. aeruginosa membranes. Additionally, we found that C10 in the water milieu can access the catalytic site of active monomeric PlaF, contributing to the competitive component of C10-mediated PlaF inhibition. Our study provides mechanistic insights into how medium-chain FFAs may regulate the activity of PlaF, a potential bacterial drug target.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Spectroscopic and QM/MM studies of the Cu(I) binding site of the plant ethylene receptor ETR1.

Herein, we present, to our knowledge, the first spectroscopic characterization of the Cu(I) active site of the plant ethylene receptor ETR1. The x-ray absorption (XAS) and extended x-ray absorption fine structure (EXAFS) spectroscopies presented here establish that ETR1 has a low-coordinate Cu(I) site. The EXAFS resolves a mixed first coordination sphere of N/O and S scatterers at distances consistent with potential histidine and cysteine residues. This finding agrees with the coordination of residues C65 and H69 to the Cu(I) site, which are critical for ethylene activity and well conserved. Furthermore, the Cu K-edge XAS and EXAFS of ETR1 exhibit spectroscopic changes upon addition of ethylene that are attributed to modifications in the Cu(I) coordination environment, suggestive of ethylene binding. Results from umbrella sampling simulations of the proposed ethylene binding helix of ETR1 at a mixed quantum mechanics/molecular mechanics level agree with the EXAFS fit distance changes upon ethylene binding, particularly in the increase of the distance between H69 and Cu(I), and yield binding energetics comparable with experimental dissociation constants. The observed changes in the copper coordination environment might be the triggering signal for the transmission of the ethylene response.

Biophysical and pharmacokinetic characterization of a small-molecule inhibitor of RUNX1/ETO tetramerization with anti-leukemic effects.

Acute myeloid leukemia (AML) is a malignant disease of immature myeloid cells and the most prevalent acute leukemia among adults. The oncogenic homo-tetrameric fusion protein RUNX1/ETO results from the chromosomal translocation t(8;21) and is found in AML patients. The nervy homology region 2 (NHR2) domain of ETO mediates tetramerization; this oligomerization is essential for oncogenic activity. Previously, we identified the first-in-class small-molecule inhibitor of NHR2 tetramer formation, 7.44, which was shown to specifically interfere with NHR2, restore gene expression down-regulated by RUNX1/ETO, inhibit the proliferation of RUNX1/ETO-depending SKNO-1 cells, and reduce the RUNX1/ETO-related tumor growth in a mouse model. However, no biophysical and structural characterization of 7.44 binding to the NHR2 domain has been reported. Likewise, the compound has not been characterized as to physicochemical, pharmacokinetic, and toxicological properties. Here, we characterize the interaction between the NHR2 domain of RUNX1/ETO and 7.44 by biophysical assays and show that 7.44 interferes with NHR2 tetramer stability and leads to an increase in the dimer population of NHR2. The affinity of 7.44 with respect to binding to NHR2 is Klig = 3.75 ± 1.22 µM. By NMR spectroscopy combined with molecular dynamics simulations, we show that 7.44 binds with both heteroaromatic moieties to NHR2 and interacts with or leads to conformational changes in the N-termini of the NHR2 tetramer. Finally, we demonstrate that 7.44 has favorable physicochemical, pharmacokinetic, and toxicological properties. Together with biochemical, cellular, and in vivo assessments, the results reveal 7.44 as a lead for further optimization towards targeted therapy of t(8;21) AML.

A MademoiseLLE domain binding platform links the key RNA transporter to endosomes.

Spatiotemporal expression can be achieved by transport and translation of mRNAs at defined subcellular sites. An emerging mechanism mediating mRNA trafficking is microtubule-dependent co-transport on shuttling endosomes. Although progress has been made in identifying various components of the endosomal mRNA transport machinery, a mechanistic understanding of how these RNA-binding proteins are connected to endosomes is still lacking. Here, we demonstrate that a flexible MademoiseLLE (MLLE) domain platform within RNA-binding protein Rrm4 of Ustilago maydis is crucial for endosomal attachment. Our structure/function analysis uncovered three MLLE domains at the C-terminus of Rrm4 with a functionally defined hierarchy. MLLE3 recognises two PAM2-like sequences of the adaptor protein Upa1 and is essential for endosomal shuttling of Rrm4. MLLE1 and MLLE2 are most likely accessory domains exhibiting a variable binding mode for interaction with currently unknown partners. Thus, endosomal attachment of the mRNA transporter is orchestrated by a sophisticated MLLE domain binding platform.

Structural, mechanistic, and physiological insights into phospholipase A-mediated membrane phospholipid degradation in Pseudomonas aeruginosa.

Cells steadily adapt their membrane glycerophospholipid (GPL) composition to changing environmental and developmental conditions. While the regulation of membrane homeostasis via GPL synthesis in bacteria has been studied in detail, the mechanisms underlying the controlled degradation of endogenous GPLs remain unknown. Thus far, the function of intracellular phospholipases A (PLAs) in GPL remodeling (Lands cycle) in bacteria is not clearly established. Here, we identified the first cytoplasmic membrane-bound phospholipase A1 (PlaF) from Pseudomonas aeruginosa, which might be involved in the Lands cycle. PlaF is an important virulence factor, as the P. aeruginosa ΔplaF mutant showed strongly attenuated virulence in Galleria mellonella and macrophages. We present a 2.0-Å-resolution crystal structure of PlaF, the first structure that reveals homodimerization of a single-pass transmembrane (TM) full-length protein. PlaF dimerization, mediated solely through the intermolecular interactions of TM and juxtamembrane regions, inhibits its activity. The dimerization site and the catalytic sites are linked by an intricate ligand-mediated interaction network, which might explain the product (fatty acid) feedback inhibition observed with the purified PlaF protein. We used molecular dynamics simulations and configurational free energy computations to suggest a model of PlaF activation through a coupled monomerization and tilting of the monomer in the membrane, which constrains the active site cavity into contact with the GPL substrates. Thus, these data show the importance of the PlaF-mediated GPL remodeling pathway for virulence and could pave the way for the development of novel therapeutics targeting PlaF.

Mapping the helix arrangement of the reconstituted ETR1 ethylene receptor transmembrane domain by EPR spectroscopy.

The plant ethylene receptor ETR1 is a key player in the perception of the phytohormone and subsequent downstream ethylene signal transmission, crucial for processes such as ripening, senescence and abscission. However, to date, there is sparse structural knowledge about the transmembrane sensor domain (TMD) of ETR1 that is responsible for the binding of the plant hormone and initiates the downstream signal transmission. Sequence information and ab initio modelling suggest that the TMD consists of three transmembrane helices. Here, we combined site-directed spin labelling with electron paramagnetic resonance spectroscopy and obtained distance restraints for liposome-reconstituted ETR1_TMD on the orientation and arrangement of the transmembrane helices. We used these data to scrutinize different computational structure predictions of the TMD.

A phospholipase B from Pseudomonas aeruginosa with activity towards endogenous phospholipids affects biofilm assembly.

Pseudomonas aeruginosa is a severe threat to immunocompromised patients due to its numerous virulence factors and biofilm-mediated multidrug resistance. It produces and secretes various toxins with hydrolytic activities including phospholipases. However, the function of intracellular phospholipases for bacterial virulence has still not been established. Here, we demonstrate that the hypothetical gene pa2927 of P. aeruginosa encodes a novel phospholipase B named PaPlaB. At reaction equilibrium, PaPlaB purified from detergent-solubilized membranes of E. coli released fatty acids (FAs) from sn-1 and sn-2 positions of phospholipids at the molar ratio of 51:49. PaPlaB in vitro hydrolyzed P. aeruginosa phospholipids reconstituted in detergent micelles and phospholipids reconstituted in vesicles. Cellular localization studies indicate that PaPlaB is a cell-bound PLA of P. aeruginosa and that it is peripherally bound to both membranes in E. coli, yet the active form was predominantly associated with the cytoplasmic membrane of E. coli. Decreasing the concentration of purified and detergent-stabilized PaPlaB leads to increased enzymatic activity, and at the same time triggers oligomer dissociation. We showed that the free FA profile, biofilm amount and architecture of the wild type and ΔplaB differ. However, it remains to be established how the PLB activity of PaPlaB is regulated by homooligomerisation and how it relates to the phenotype of the P. aeruginosa ΔplaB. This novel putative virulence factor contributes to our understanding of phospholipid degrading enzymes and might provide a target for new therapeutics against P. aeruginosa biofilms.

Unsaturated fatty acids augment protein transport via the SecA:SecYEG translocon.

The translocon SecYEG and the associated ATPase SecA form the primary protein secretion system in the cytoplasmic membrane of bacteria. The secretion is essentially dependent on the surrounding lipids, but the mechanistic understanding of their role in SecA : SecYEG activity is sparse. Here, we reveal that the unsaturated fatty acids (UFAs) of the membrane phospholipids, including tetraoleoyl-cardiolipin, stimulate SecA : SecYEG-mediated protein translocation up to ten-fold. Biophysical analysis and molecular dynamics simulations show that UFAs increase the area per lipid and cause loose packing of lipid head groups, where the N-terminal amphipathic helix of SecA docks. While UFAs do not affect the translocon folding, they promote SecA binding to the membrane, and the effect is enhanced up to fivefold at elevated ionic strength. Tight SecA : lipid interactions convert into the augmented translocation. Our results identify the fatty acid structure as a notable factor in SecA : SecYEG activity, which may be crucial for protein secretion in bacteria, which actively change their membrane composition in response to their habitat.

Substrate Access Mechanism in a Novel Membrane-Bound Phospholipase A of Pseudomonas aeruginosa Concordant with Specificity and Regioselectivity.

PlaF is a cytoplasmic membrane-bound phospholipase A1 from Pseudomonas aeruginosa that alters the membrane glycerophospholipid (GPL) composition and fosters the virulence of this human pathogen. PlaF activity is regulated by a dimer-to-monomer transition followed by tilting of the monomer in the membrane. However, how substrates reach the active site and how the characteristics of the active site tunnels determine the activity, specificity, and regioselectivity of PlaF for natural GPL substrates have remained elusive. Here, we combined unbiased and biased all-atom molecular dynamics (MD) simulations and configurational free-energy computations to identify access pathways of GPL substrates to the catalytic center of PlaF. Our results map out a distinct tunnel through which substrates access the catalytic center. PlaF variants with bulky tryptophan residues in this tunnel revealed decreased catalysis rates due to tunnel blockage. The MD simulations suggest that GPLs preferably enter the active site with the sn-1 acyl chain first, which agrees with the experimentally demonstrated PLA1 activity of PlaF. We propose that the acyl chain-length specificity of PlaF is determined by the structural features of the access tunnel, which results in favorable free energy of binding of medium-chain GPLs. The suggested egress route conveys fatty acid (FA) products to the dimerization interface and, thus, contributes to understanding the product feedback regulation of PlaF by FA-triggered dimerization. These findings open up opportunities for developing potential PlaF inhibitors, which may act as antibiotics against P. aeruginosa.

TopProperty: Robust Metaprediction of Transmembrane and Globular Protein Features Using Deep Neural Networks.

Transmembrane proteins (TMPs) are critical components of cellular life. However, due to experimental challenges, the number of experimentally resolved TMP structures is severely underrepresented in databases compared to their cellular abundance. Prediction of (per-residue) features such as transmembrane topology, membrane exposure, secondary structure, and solvent accessibility can be a useful starting point for experimental design or protein structure prediction but often requires different computational tools for different features or types of proteins. We present TopProperty, a metapredictor that predicts all of these features for TMPs or globular proteins. TopProperty is trained on datasets without bias toward a high number of sequence homologs, and the predictions are significantly better than the evaluated state-of-the-art primary predictors on all quality metrics. TopProperty eliminates the need for protein type- or feature-tailored tools, specifically for TMPs. TopProperty is freely available as a web server and standalone at https://cpclab.uni-duesseldorf.de/topsuite/.

Computational Analyses of the AtTPC1 (Arabidopsis Two-Pore Channel 1) Permeation Pathway.

Two Pore Channels (TPCs) are cation-selective voltage- and ligand-gated ion channels in membranes of intracellular organelles of eukaryotic cells. In plants, the TPC1 subtype forms the slowly activating vacuolar (SV) channel, the most dominant ion channel in the vacuolar membrane. Controversial reports about the permeability properties of plant SV channels fueled speculations about the physiological roles of this channel type. TPC1 is thought to have high Ca2+ permeability, a conclusion derived from relative permeability analyses using the Goldman-Hodgkin-Katz (GHK) equation. Here, we investigated in computational analyses the properties of the permeation pathway of TPC1 from Arabidopsis thaliana. Using the crystal structure of AtTPC1, protein modeling, molecular dynamics (MD) simulations, and free energy calculations, we identified a free energy minimum for Ca2+, but not for K+, at the luminal side next to the selectivity filter. Residues D269 and E637 coordinate in particular Ca2+ as demonstrated in in silico mutagenesis experiments. Such a Ca2+-specific coordination site in the pore explains contradicting data for the relative Ca2+/K+ permeability and strongly suggests that the Ca2+ permeability of SV channels is largely overestimated from relative permeability analyses. This conclusion was further supported by in silico electrophysiological studies showing a remarkable permeation of K+ but not Ca2+ through the open channel.

Interdependence of a mechanosensitive anion channel and glutamate receptors in distal wound signaling.

Glutamate has dual roles in metabolism and signaling; thus, signaling functions must be isolatable and distinct from metabolic fluctuations, as seen in low-glutamate domains at synapses. In plants, wounding triggers electrical and calcium (Ca2+) signaling, which involve homologs of mammalian glutamate receptors. The hydraulic dispersal and squeeze-cell hypotheses implicate pressure as a key component of systemic signaling. Here, we identify the stretch-activated anion channel MSL10 as necessary for proper wound-induced electrical and Ca2+ signaling. Wound gene induction, genetics, and Ca2+ imaging indicate that MSL10 acts in the same pathway as the glutamate receptor­like proteins (GLRs). Analogous to mammalian NMDA glutamate receptors, GLRs may serve as coincidence detectors gated by the combined requirement for ligand binding and membrane depolarization, here mediated by stretch activation of MSL10. This study provides a molecular genetic basis for a role of mechanical signal perception and the transmission of long-distance electrical and Ca2+ signals in plants.

The Puzzle of Metabolite Exchange and Identification of Putative Octotrico Peptide Repeat Expression Regulators in the Nascent Photosynthetic Organelles of Paulinella chromatophora.

The endosymbiotic acquisition of mitochondria and plastids more than one billion years ago was central for the evolution of eukaryotic life. However, owing to their ancient origin, these organelles provide only limited insights into the initial stages of organellogenesis. The cercozoan amoeba Paulinella chromatophora contains photosynthetic organelles-termed chromatophores-that evolved from a cyanobacterium ∼100 million years ago, independently from plastids in plants and algae. Despite the more recent origin of the chromatophore, it shows tight integration into the host cell. It imports hundreds of nucleus-encoded proteins, and diverse metabolites are continuously exchanged across the two chromatophore envelope membranes. However, the limited set of chromatophore-encoded solute transporters appears insufficient for supporting metabolic connectivity or protein import. Furthermore, chromatophore-localized biosynthetic pathways as well as multiprotein complexes include proteins of dual genetic origin, suggesting that mechanisms evolved that coordinate gene expression levels between chromatophore and nucleus. These findings imply that similar to the situation in mitochondria and plastids, also in P. chromatophora nuclear factors evolved that control metabolite exchange and gene expression in the chromatophore. Here we show by mass spectrometric analyses of enriched insoluble protein fractions that, unexpectedly, nucleus-encoded transporters are not inserted into the chromatophore inner envelope membrane. Thus, despite the apparent maintenance of its barrier function, canonical metabolite transporters are missing in this membrane. Instead we identified several expanded groups of short chromatophore-targeted orphan proteins. Members of one of these groups are characterized by a single transmembrane helix, and others contain amphipathic helices. We hypothesize that these proteins are involved in modulating membrane permeability. Thus, the mechanism generating metabolic connectivity of the chromatophore fundamentally differs from the one for mitochondria and plastids, but likely rather resembles the poorly understood mechanism in various bacterial endosymbionts in plants and insects. Furthermore, our mass spectrometric analysis revealed an expanded family of chromatophore-targeted helical repeat proteins. These proteins show similar domain architectures as known organelle-targeted expression regulators of the octotrico peptide repeat type in algae and plants. Apparently these chromatophore-targeted proteins evolved convergently to plastid-targeted expression regulators and are likely involved in gene expression control in the chromatophore.

Structural Model of the ETR1 Ethylene Receptor Transmembrane Sensor Domain.

The structure, mechanism of action and copper stoichiometry of the transmembrane sensor domain of the plant ethylene receptor ETR1 and homologs have remained elusive, hampering the understanding on how the perception of the plant hormone ethylene is transformed into a downstream signal. We generated the first structural model of the transmembrane sensor domain of ETR1 by integrating ab initio structure prediction and coevolutionary information. To refine and independently validate the model, we determined protein-related copper stoichiometries on purified receptor preparations and explored the helix arrangement by tryptophan scanning mutagenesis. All-atom molecular dynamics simulations of the dimeric model reveal how ethylene can bind proximal to the copper ions in the receptor, illustrating the initial stages of the ethylene perception process.

Calcium-Promoted Interaction between the C2-Domain Protein EHB1 and Metal Transporter IRT1 Inhibits Arabidopsis Iron Acquisition.

Iron is a key transition element in the biosphere and is crucial for living organisms, although its cellular excess can be deleterious. Maintaining the balance of optimal iron availability in the model plant Arabidopsis (Arabidopsis thaliana) requires the precise operation of iron import through the principal iron transporter IRON-REGULATED TRANSPORTER1 (IRT1). Targeted inhibition of IRT1 can prevent oxidative stress, thus promoting plant survival. Here, we report the identification of an IRT1 inhibitor, namely the C2 domain-containing peripheral membrane protein ENHANCED BENDING1 (EHB1). EHB1 interacts with the cytoplasmically exposed variable region of IRT1, and we demonstrate that this interaction is greatly promoted by the presence of calcium. We found that EHB1 binds lipids characteristic of the plasma membrane, and the interaction between EHB1 and plant membranes is calcium-dependent. Molecular simulations showed that EHB1 membrane binding is a two-step process that precedes the interaction between EHB1 and IRT1. Genetic and physiological analyses indicated that EHB1 acts as a negative regulator of iron acquisition. The presence of EHB1 prevented the IRT1-mediated complementation of iron-deficient fet3fet4 yeast (Saccharomyces cerevisiae). Our data suggest that EHB1 acts as a direct inhibitor of IRT1-mediated iron import into the cell. These findings represent a major step in understanding plant iron acquisition, a process that underlies the primary production of bioavailable iron for land ecosystems.

PACKMOL-Memgen: A Simple-To-Use, Generalized Workflow for Membrane-Protein-Lipid-Bilayer System Building.

We present PACKMOL-Memgen, a simple-to-use, generalized workflow for automated building of membrane-protein-lipid-bilayer systems based on open-source tools including Packmol, memembed, pdbremix, and AmberTools. Compared with web-interface-based related tools, PACKMOL-Memgen allows setup of multiple configurations of a system in a user-friendly and efficient manner within minutes. The generated systems are well-packed and thus well-suited as starting configurations in MD simulations under periodic boundary conditions, requiring only moderate equilibration times. PACKMOL-Memgen is distributed with AmberTools and runs on most computing platforms, and its output can also be used for CHARMM or adapted to other molecular-simulation packages.

Nutrient exchange in arbuscular mycorrhizal symbiosis from a thermodynamic point of view.

To obtain insights into the dynamics of nutrient exchange in arbuscular mycorrhizal (AM) symbiosis, we modelled mathematically the two-membrane system at the plant-fungus interface and simulated its dynamics. In computational cell biology experiments, the full range of nutrient transport pathways was tested for their ability to exchange phosphorus (P)/carbon (C)/nitrogen (N) sources. As a result, we obtained a thermodynamically justified, independent and comprehensive model of the dynamics of the nutrient exchange at the plant-fungus contact zone. The predicted optimal transporter network coincides with the transporter set independently confirmed in wet-laboratory experiments previously, indicating that all essential transporter types have been discovered. The thermodynamic analyses suggest that phosphate is released from the fungus via proton-coupled phosphate transporters rather than anion channels. Optimal transport pathways, such as cation channels or proton-coupled symporters, shuttle nutrients together with a positive charge across the membranes. Only in exceptional cases does electroneutral transport via diffusion facilitators appear to be plausible. The thermodynamic models presented here can be generalized and adapted to other forms of mycorrhiza and open the door for future studies combining wet-laboratory experiments with computational simulations to obtain a deeper understanding of the investigated phenomena.

The Nature of the Donor Motif in Acceptor-Bridge-Donor Dyes as an Influence in the Electron Photo-Injection Mechanism in DSSCs.

The combination and balance of acceptor(A)-bridge-donor(D) architecture of molecules confer suitable attributes and/or properties to act as efficient light-harvesting and sensitizers in dye sensitized solar cells (DSSCs). An important process in a DSSC performance is the electron photoinjection (PI) mechanism which can take place either via type I (indirect), that consists in injecting from the excited state of the dye to the semiconductor, or type II (direct), where the PI is from the ground state of the dye to the semiconductor upon photoexcitation. Here, we present a computational study about the role of the donor motif in the PI mechanisms displayed from a family of 11 A-bridge-D structured dyes to a (TiO2)15 anatase cluster. To this end, different donor motifs (D1-D11) were evaluated while the A and bridge motifs remained the same. All the computations were carried out within the DFT framework, using the B3LYP, PW91, PBE, M06L and CAM-B3LYP functionals. The 6-31G(d) basis set was employed for nonmetallic atoms and the LANL2DZ pseudopotential for Ti atoms. The solvation effects were incorporated using the polarized continuum model (PCM) for acetonitrile. As benchmark systems, alizarin and naphthalenediol dyes were analyzed, as they are known to undergo Type I and Type II PI pathways in DSSCs, respectively. Donors in the studied family of dyes could influence to drive Type I or II PI since it was found that D2 could show some Type II PI route, showing a new absorption band, although with CAM-B3LYP this shows a very low oscillator strength, while the remaining dyes behave according to Type I photoinjectors. Finally, the photovoltaic parameters that govern the light absorption process were evaluated, as the use of these criteria could be applied to predict the efficiency of the studied dyes in DSSCs devices.